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ABSTRACT 


A computer technique is proposed for a simple practical method of auto- 
matically designing tower structures. Dynamic programming is used to find 
the optimum geometric configuration of the structural members , while the 
member sizes are proportioned by direct iteration. 

Tower structures are particularly suited to this method of automatic 
design since the rapidity of the analysis and design depends primarily upon 
substructuring. Substructuring of towers is comparatively simple because 
interaction between adjacent substructures can be simulated with reasonable 
accuracy. Typical examples are presented to illustrate the method. 


INTRODUCTION 


Dynamic programming has long been recognised as an extremely powerful 
optimization technique, particularly for problems of a discontinuous nature. 
High-dimensional problems, however, result in a large amount of computation, 
but this can be reduced by a successive approximation method. 

A 3-dimensional Dynamic Programming Successive Approximation technique 
is used here to obtain an optimum (least weight) geometrical configuration 
for the design of tower structures. Concurrently, a simple direct itera- 
tion procedure is used to select the optimum member sizes from any list of 
section properties provided. 

If the structure were to be designed as a whole, a change of geometric 
configuration would need a re-analysis. This means that a large amount 
of computation would be required. 

Substructuring has been introduced to reduce the overall problem into 
smaller stages so that the analysis can be performed more rapidly. 

If any sections recommended by the computer are not desirable for 
practical reasons the re-input of a complete new set of data is unnecessary. 
Only the list(s) of sections and the data relating to section types need 
be changed. 


L 
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Three examples of tower structures are given in this paper: 

1. A sample plane-truss tower 

2. A comparison of the optimum weight of a triangular tower referred to by 
Kuzmanovic, Willems and Thomas (ref. l) 

3. A typical transmission tower used in South Africa 

DYNAMIC PROGRAMMING - BASIC CONCEPTS 


The computational effort required to find all the possible solutions for 
large problems can become unmanageable without the use of Dynamic Programming. 
This technique bypasses this problem by considering the possible decisions 
to be made at each stage of the solution. 

Dynamic Programming can be described as a technique for methodically 
selecting an optimum solution of a multi-stage decision problem. This 
mathematical technique can be used for problems where a sequence of decisions 
are dependent upon one another and each decision influences the system's 
response to future decisions. A set of solutions can be categorized so that 
one may be judged to be better than another in some pre-defined manner. 

In a sequence of decisions, the current state of the sequence is assessed 
as f(x(k-l)) and the succeeding one as f (x(k) ) . Without considering the 

whole chain of past and future decisions, except that they contribute to 
f(x(k-l)), the best decision can be found from: 

f(x(k)) = min[t(x(k); x(k-l)) + f(x(k-l))] 

where t is the assessment between stages k-1 and k and 
x is a state variable 

Dynamic Programming is based on a repeated use of this idea. 

A simple network problem posed by Bellman and Dreyfus (refs. 2 & 3) and 
discussed by Palmer (ref. 9) demonstrates the process excellently. In the 
course of the solution of this problem, two central ideas are used. The 
first is the idea of imbedding; this means that the overall optimization 
problem consists of a number of smaller problems imbedded within the whole, 
which can be solved independently. In the second idea, an optimum solution 
can be found from a sequence of decisions by imagining that the final solution 
is broken up into a series of simpler decisions. 

Problems of this type become very tedious when the order of the state 
variables (x(k)) becomes large. For example, various state variable vectors 
(denoted x (k)) may be present at any particular stage k of the calculation, 
which compSicates the matter. For this reason, the dimension (n) of the 
problem is decreased to a single variable vector by the use of the Dynamic 
Programming Successive Approximation (DPSA) technique. 
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The geometric design of a tower structure can he degenerated into a 
series of simple decision processes as prescribed by the DPSA technique by 
using substructuring. 


SUBSTRU C TURING 


The benefit of substructuring in the design of towers is that the optimum 
section sizes can be determined rapidly within each substructure. The inter- 
action between adjacent substructures is comparatively simple and can be 
simulated with reasonable ease which makes this technique particularly 
attractive . 

The geometry at the interface of each substructure is uniquely defined 
by the coordinates of the member ends which are 'cut' at the interfaces. 

Hence, in Figure 1, the coordinates at the right hand side of substructure 1 

are defined by x q , y ; x^, y^ and for substructure 2 by x^, y^; x^, y^ 

and so on for the other substructures. 

The assumption is made that the forces transferred from one substructure 

to the next can be calculated from statics. For example, the interaction 

forces at the interface between substructures 1 and 2 are found by applying 

equivalent loads and moments at the central point A; i.e., the vertical 

force P^ is resolved into a load P and a moment P x h, while the horizontal 

force P is resolved into a load P. and a moment P, x (y« - v, ) . 

h h h J 4 J 1 

Each pin-jointed substructure is analysed by the displacement method, 
which requires the equivalent loads and moments at A to be applied as point 
loads at the interface nodes (x^, y^) and ( -x^ , y ). The substructures are 
now analysed and designed independently assuming that the lower interface 
nodes are constrained. The member sizes are selected by a simple direct 
iteration procedure. 


DIRECT ITERATION 


A list of sections is provided so that the actual member sizes required 
in each substructure can be selected automatically by the direct iteration 
procedure. The calculated stresses are compared with permissible stresses 
so that the ratio between them is a minimum. The substructure is reanalysed 
each time the member sizes are revised until the results from successive 
iterations are identical. 


THE DYNAMIC PROGRAMMING SUCCESSIVE APPROXIMATION (DPSA) TECHNIQUE 


The coordinates of the interface nodes in Figure 1 can be changed at any 
stage and the members designed accordingly by the direct iteration method. 
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Consequently, the weight of each configuration of members in all substructures 
can be calculated. The optimum solution is then the combination of possible 
configurations which results in the least weight of the total structure. The 
dynamic programming technique sets about this calculation in an organized 
methodical way. 

Let us consider a tower similar to that in Figure 1. The interface 
nodes, which define the configuration of the tower, can be positioned ‘in space 
by a set of three-dimensional coordinate values. Let the structure be 
symmetrical about the XZ and YZ planes; therefore a node placed in the first 
octant will define the shape of the tower at that level. Let the coordinates 
of the interface node be (x (k) , Xg(k), x^(k)). This corresponds to the 
x, y and z coordinates of the primary node of substructure 1, where k also 
denotes the upper interface level of that substructure and k = 0 represents 
the ground level. The n-dimensional DPSA method is therefore well suited 
to structures of this type, where n = 3. For example, at level k = 2 some 
possible values of the state variables are shown in Table 1. The control 
variables u^Ck) shown correspond to the identification number in each set. 

The values of u(k) (i.e. u.,(k), u_(k), u^(k)) and hence the values of 

x (k) , x (k), x^(k) at each level k must be found so that the overall weight 
or the structure is a minimum. 

The DPSA method requires an initial solution to the problem. The average 
value of the state variables at each interface is a suitable initial solution. 
To proceed, only one state variable is altered, while the remainder stay at 
their initial values. In this way a single-dimensional dynamic programming 
procedure with respect to this variable is carried out. The process is then 
continued, the new value of the first variable is retained and one of the other 
state variables is altered. All the variables are processed in this way; 
the first cycle is complete when all the state variables have been altered. 
Subsequent cycles follow the same procedure as above. The DPSA method has 
converged to its optimum solution when no weight difference is recorded 
between successive cycles. This method is found to converge rapidly to an 
optimum weight solution. The structure comprises a number of substructures. 
The geometric configuration between interfaces is regarded as a substructure. 

The direct iteration method is used to find the most satisfactory 
structural design and hence the weight for every geometric configuration of 
each substructure. A least weight path is followed through all the sub- 
structures to determine the optimum configuration of the total structure. The 
explanation of the DPSA method can be simplified by the following elementary 
example . 

Example : Consider a simple 2-dimensional tower which consists of 2 

substructures - Figure 2. 

Let the state variables x(l) and x(2) vary in three steps on each side 
of the vertical axis of symmetry and the other state variables i.e. y(k) be 
kept constant. In addition, let the state variable x(k) at k = 0 be 
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constrained to a single value. The accumulated cost (or weight) of the sub- 
structures up to level k for each position of x can he denoted by: 

Ch(k) 

where i denotes the position of x(k) on interface k. 

Let the cost of a single substructure k be denoted by: 

t 13 (k) 

where i and j denote the positions of the state variables x(k) for sub- 
structure k at the upper and lower interfaces respectively. 

Figure 3 shows all the possible geometrical configurations within the 
constraints given. 

The costs of each of the 3 configurations in the first substructure, due 
to the varying values of x(l), are calculated by: 

C-^l) = t u (l) 

c 2 d) = t 21 (l) 

C 3 (l) = t 31 d) 
and are shown in Figure U . 

Similarly the values of t. .(2) can be obtained for the second sub- 
structure. The least accumulated weight at point i on interface k can be 
found from: 

C. (2) = min[t ..(2) + C.(l)] for j = 1 to 3 

^-0 J 

The optimum configuration of the structure at level 2 is determined hy choos- 
ing the least value of C.(2). The optimum path can then be traced back 
through the calculations to find the optimum configuration of the entire 
structure . 

A computer program, which uses the ideas discussed above, has been written 
to design general tower structures. Three examples included are: 
a. A simple 3-cell plane-truss tower. b. A 3-legged transmission tower, 
c. A practical, rectangular plan transmission tower. 

Example 1. 3-substructure plane-truss tower 

The tower shown in Figure 5 supports two loads at the upper level of 
-15 kN in the y-direction and -10 kN in the x-direction at nodes 7 and 8 
respectively. The possible dimensions at the k levels are, in x-direction: 



1 

2 

3 

u 

5 

Base 

1,6 

1,8 

2,0 

2,2 

2,k 

Level 1 

1,1 

1,3 

1,5 

1,7 

1,9 

Level 2 

0,8 

0,9 

1,0 

1,1 

1,2 

Level 3 

0,5 

- 

- 

- 

- 


and in z-direction: 
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1 2 

3 

4 5 


Base 

0,0 

- 

— — 


Level 1 

1,8 1,9 

2,0 

2,1 2,2 


Level 2 

3,8 3,9 

4,0 

4,1 4,2 


Level 3 

6,0 

- 



Each substructure is to be 

designed using the following 

sections : 

1. Upright member - channel sections. 

2. Diagonal member 

- pipe sections 

3. Horizontal 

member - angle sections. 



Results : 

The final shape 

is shown 

in Figure 6 . 


Structural 

Design : 




Substructure 

Members 


Size - mm 

Type 

1 

a ,d 


100 x 50 x 6 

Channel 


b ,c 


38 ,l<j> x 2 

Pipe 


e 


40 x 40 x 3 

Angle 

2 

f,i 


100 x 50 x 6 

Channel 


g,h 


38 , l(f> x 2 

Pipe 


j 


25 x 25 x 3 

Angle 

3 

k,n 


76 x 38 x 5,1 

Channel 


i ,m 


48,54) x 2 

Pipe 


0 


25 x 25 x 3 

Angle 

Total weight of 

structure, l,46l 

kN; total computer time, 28, 

39 sec UNIVAC 1106 


Example 2. A three-legged transmission tower shown in Figure 7 is 
considered. It contains 56 nodes and 175 members. The loads used correspond 
to those used by Kuzmanovic et al. (ref. l). 

Load case 1: Basic wind free in transverse direction 


Position Load kN Direction 

A,B 26,7 - z 

8,9 - x 

C 80,0 - z 

31.0 - x 

D,E 15,6 - x 

40.0 - z 

Load case 2: 0,707 Basic wind in the transverse direction 

0,707 Basic wind in the longitudinal direction 

Position Load kN Direction 

A,B l6,0 - z 

5,4 - x 

C 48,0 - z 

18,7 - x 

D,E 9,4 - x 

11,6 + y 

24.0 - z 
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Load case 3: No wind 



Position 

Load kN 

Direction 

A,B 

32,0 

- z 

C 

160,0 

- z 

D,E 

80,0 

- z 


Member Groups : Table 2 shows the member groups used in each substructure . 

Angle sizes: Table 3 shows the list of angles which were used to 

produce a structural design. 


Results: The structural design for each load case is shown in Figure 8. 
The horizontal axis represents, in ascending order, -the sizes of angle avail- 
able for the design. The vertical scale represents the member groups within 
each substructure. Computer times and structure weights are: 


Load case 


1 

2 

3 


No. of state variable 
positions at each level 
Horizontal Vertical 

(radial) (elevation) 

5 1 

5 1 

5 1 


Weight kN Computer 

time 

UNIVAC 1106 


23,256 6 min 27 sec 

16,786 4 min 12,17 sec 

17,428 3 min 36,89 sec 


The optimum weight found by Kuzmanovic et al. (ref. l) was 23,51 kN. The 
weights above compare favourably with this value. 


Example 3 . A practical transmission tower similar to those currently 
in use in South Africa is considered. It contains 115 nodes and 336 members. 
The shape of tower is shown in Figure 9. An equivalent set of loads have 
been calculated from those used in Example 2. Angle sizes used in the design 
are shown in Table 3 . 

Member Groups: Table 4 shows the member groups used in each substructure. 


Results: The structural design for each load case is shown in Figure 10. 


Computer times and structural weights are: 

No. of state variable 
positions at each level 


Load case 


1 

2 


x-dir. y-dir. 


3 

3 


z-dir. 

1 

1 


Weight kN 


57,877 

48 , 541 * 


Computer 

time 

UNIVAC 1106 

46 min 04,252 sec 
31 min 15,069 sec 


CONCLUSION 


This method of automatically determining the optimum geometric configura- 
tion and member sizes substantially reduces the effort involved in the design 
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of tower structures . The dynamic programming successive approximations 
technique is effective for structures of this type, where substructuring pro- 
vides the necessary static variables. The displacement method is admirably 
suited to the analysis of these structures. Direct iteration is the simplest 
method of selecting the optimum member sizes from any given list of section 
properties. Three examples have been discussed. In the 15 member plane- 
truss example considered as 3 substructures, the computational time taken on 
a UNIVAC 1106 is 28,39 seconds. A waisted tower shape of structure is the 
optimum geometric solution. Example 2 is a three-legged tower with 175 
members and 56 nodes. The computational time with two state variables is 
6 minutes 27 seconds. The weight compares favourably with that given in 
reference 11 by Kuzmanovic et al. Example 3 is a practical transmission tower 
similar to those currently used in South Africa. This structure consists 
of 336 members and 115 nodes and the computational time required for an 
optimum feasible solution is b6 minutes per load case. 
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TABLE 1 


A POSSIBLE SET OF STATE AND CONTROL VARIABLE VECTORS 


STATE 

VARIABLES 


x x (2) 

x 2 (2) 

x 3 (2) 

0,8 

0,7 

0,9 

0,9 

0,75 

0,95 

1,0 

0;8 

1,0 

1,1 

0,85 

1,05 

1,2 

0,9 

1,1 


IDENTIFICATION NO u (k) 
n = 1 to 3 n 
u(2) 

1 

2 

3 

h 

5 


TABLE 2 - 

Type 

Main leg members , a 
Horizontal leg members, b 
Diagonal leg members, c 
Main leg members, d,g,j,m,p 
Horizontal members, e,h,k,n,q 
Diagonal members, f,i,£,o,r 


MEMBER GROUPS 

Substructures 
12 3 4 



5 


1 

2 

3 


6 


1 

2 

3 
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TABLE 3 


AVAILABLE ANGLE SIZES 


Number 

Size ■ 

- ram 

Number 


Size ■ 

- mm 

1 

25 

X 

25 

X 

3 

11 

80 

X 

80 

x 8 

2 

30 

X 

30 

X 

3 

12 

90 

X 

90 

x 8 

3 

4o 

X 

4o 

X 

3 

13 

100 

X 

100 

x 8 

4 

45 

X 

45 

X 

3 

i4 

100 

X 

100 

x 10 

5 

45 

X 

45 

X 

5 

15 

120 

X 

120 

x 10 

6 

50 

X 

50 

X 

5 

16 

120 

X 

120 

CVI 

H 

W 

7 

60 

X 

60 

X 

5 

17 

150 

X 

150 

x 12 

8 

60 

X 

60 

X 

6 

18 

150 

X 

H 

\s\ 

0 

H 

CO 

9 

70 

X 

70 

X 

6 

19 

200 

X 

200 

X 16 

10 

80 

X 

80 

X 

6 

20 

200 

X 

200 

X 24 


TABLE 4 - MEMBER GROUPS 


Type 

Main leg members, a,f,j 
Secondary leg members, b,g,k 
Leg horizontal members, c 
Diagonal members, d 
Interface members, e,i,m 
Leg bracing members, h,S, 

Main arm members, n 
Secondary arm members , o 
Bracing, p,q,r,s,t,u 
Boom members, v,w 


Substructures 

12 3 4 

111 
2 2 2 

3 

4 

5 4 4 

3 3 

1 

2 

3, 4, 5, 6, 7,8 

9,10 


111 




Figure 1.- Tower structure with substructures. 



Figure 2.- Tower structure (example). 


POSITION 



Figure 3.- Possible configurations. 


POSITION 



Figure 4.- Three cost values. 








SUBSTRUCTURE 

NUMBER 

3 


2 


1 


MEMBER 

GROUP 

NUMBER 



Figure 8.- Example 2: Structural design. 




Figure 9.- Typical tower. 
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